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Abstract. We consider the massive wave equation on asymptotically AdS spaces. We 
show that the timelike X behaves like a finite timelike boundary, on which one may impose 
the equivalent of Dirichlet, Neumann or Robin conditions for a range of (negative) mass 
parameter which includes the conformally coupled case. We demonstrate well posedness 
for the associated initial-boundary value problems at the level of regularity. We also 
prove that higher regularity may be obtained, together with an asymptotic expansion 
for the field near X. The proofs rely on energy methods, tailored to the modified energy 
introduced by Breitenlohner and Freedman. We do not assume the spacetime is stationary, 
nor that the wave equation separates. 



Among the solutions of Einstein's general theory of relativity, the maximally symmetric 
spacetimes hold a privileged position. Owing to the high level of symmetry, they serve as 
plausible 'ground states' for the gravitational field, so there is great interest in spacetimes 
which approach a maximally symmetric spacetime in some asymptotic region. Such a 
spacetime would represent an 'isolated gravitating system'. Historically, asymptotically 
flat spacetimes have been the most studied, however, recently there has been great interest 
in the asymptotically anti-de Sitter (AdS) spacetimes motivated by the putative AdS/CFT 
correspondence [1] . In the study of classical General Relativity, there have also been some 
very interesting recent results regarding the question of black hole stability [2 [3l |3] for 
asymptotically AdS black holes. 

The asymptotically AdS spacetimes approach (the covering space of) the spacetime of 
constant sectional curvature — p-, which we shall refer to as the AdS spacetime. In so called 
'global coordinates', the metric takes the form 



In contrast to the Minkowski spacetime, AdS has a timelike conformal boundary, X. Ac- 
cordingly one expects that in order to have a well posed time evolution for the equations 
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of physics in this background it is necessary to specify some boundary condition on X. In 
pioneering work [5], Breitenlohner and Freedman considered the massive wave equation 

(1.1) Dg^-^Cl) = 0, 

on a fixed anti-de Sitter background of constant sectional curvature —p- They were able 
to solve the wave equation by separation of variables, making use of the 50(2, 3) symmetry 
of AdS. The second order ordinary differential equation governing the radial part of the 
wave equation has a regular singular point at infinity, hence the field has an expansion near 
infinity: 



where X± = | it y^l + A. When the mass parameter is in the range — | < A < 0, both 
branches decay towards infinity. For a well posed problem it is necessary to place some 
constraints on the functions tp^. The usual choice would be to insist that ip' = 0, which is 
analogous to imposing a Dirichlet condition at I. This corresponds to requiring a solution 
of finite energy (we shall elaborate on this point later). Breitenlohner and Freedman showed 
that for — I < A < — I the wave equation can also be solved on the exact anti-de Sitter 
space under the assumption that ip'^ = 0, analogous to a Neumann condition. From now 
on, we will assume A to lie in this range. Breitenlohner and Freedman also introduced a 
modified or renormalised energy which is finite for both branches of the solution. As in 
the case of a finite domain, the Dirichlet and Neumann boundary conditions are not the 
only possible choices. We summarise some other possible boundary conditions in Table [l] 
below (the list is not exhaustive). 



(1) 


Dirichlet 


^p- =0 


(2) 


Inhomogeneous Dirichlet 




(3) 


Neumann 


V'+ = 


(4) 


Inhomogeneous Neumann 


V;+ = / 


(5) 


Robin 


^+ + = 


(6) 


Inhomogeneous Robin 





Table 1. Possible boundary conditions. /, (3 are functions on X. 



+ 



i^-{e,(t>,t) + o 



The homogeneous conditions (1), (3), and (5) were considered by Ishibashi and Wald 
[6l[7], who showed that they give rise to a well defined unitary evolution for the scalar wave, 
Maxwell and gravitational perturbation equations in the exact anti-de Sitter spacetime. 
This work has been extended to the Dirac equation in the work of Bachelot [8]. These 
papers use methods based on self-adjoint extensions of the elliptic part of the wave operator. 
They make crucial use of properties of the exact AdS space (in particular staticity and 
separability) which are not shared by the general class of asymptotically AdS spaces. For 
Dirichlet boundary conditions (1), (2), the work of Holzegel [9j and Vasy [lOj provides 
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well posedness results for a more general class of asymptotically AdS spaces for the range 
A > — I based on energy space methods, but does not treat the other possible boundary 
conditions. 

The aim of this work is to treat all of the boundary conditions (l)-(6), without making 
any assumptions regarding the staticity or separability of the metric. The natural ap- 
proach is that of energy estimates, however, we are forced to confront the problem that for 
boundary conditions (2)-(6) we cannot expect the standard energy to be finite. To deal 
with this, we use the renormalised energy of [5j- We show that the natural Hilbert spaces 
associated with this energy are generalisations of the usual Sobolev spaces where 'twisted' 
derivatives of the form 

OXi 

are supposed to exist in a distributional sense and belong to an appropriate space, for 
some appropriately chosen /j, related to the distance to the boundary. 

To simplify the analysis at the expense of losing the geometrical structure, we map 
the problem to a more general problem in a finite region of M" which is very closely 
analogous to that of a finite initial boundary value problem (IB VP). We introduce weak 
formulations at the level for all of the boundary conditions, and are able to show 
that these weak problems admit a unique solution. The method used to show existence 
is to approximate the problem by a suitable hyperbolic IBVP in a finite cylinder, and let 
the cylinder approach T. We then use energy estimates to extract a weakly convergent 
subsequence, in this way, we can show well posedness results for (l)-(6). We are also able 
to recover higher regularity for the solution, if more assumptions are made on the data. 
We further provide an asymptotic expansion for the solutions near infinity. 

The paper will be structured as follows. We first define the asymptotically AdS spaces 
we consider in ^ We then introduce the modified energy in f|3]and use it in ^to motivate 
weak formulations of the Dirichlet and Neumann problems, which we then show to be well 
posed. In ^we show that under stronger assumptions on data improved regularity can be 
obtained, together with the asymptotic behaviour of the solution. Finally, in f|6]we discuss 
briefiy the inhomogeneous and Robin boundary conditions and remark on the connection 
to methods involving self-adjoint extensions. We assume throughout a degree of familiarity 
with the theory of the finite IBVP, as developed for example in |TH I12j. 

Acknowledgements. I would like to thank Gustav Holzegel for introducing me to this 
problem, and for helpful comments. I would also like to thank Mihalis Dafermos, Ju- 
lian Sonner and Pau Figueras for comments. I would like to acknowledge funding from 
PIMS and NSERC. The early stages of this project were supported by Queens' College, 
Cambridge. 

2. Asymptotically AdS spaces 

We take our definition of the asymptotically AdS spaces from [9], with some slight 
modifications. In particular, we find it convenient to assume that the inner boundary is 
null. 
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Figure 1 . The asymptotic region T> 



Definition 1. We say that the smooth Lorentzian manifold {T>,g) is an asymptotically 
AdS spacetime patch with cosmological constant A = —p if the following holds 

(i) V is topologically [0, T] x [i?, oo) x S"^ and is covered by a coordinate system (r, t, x"^, y*), 
where x*, (i = 1,2 are two sets of coordinates which together cover the sphere). The 
hypersurfaces of constant t foliate T) and are everywhere spacelike, we denote them 
by Tif Each is foliated by 2— spheres S^^. The boundary of the region T> is given 
by Sq (in the past), Y^t (in the future) and finally the null hypersurface Bq generated 
by evolving Sq along a set of null geodesies directed away from the boundary such 
that any timelike geodesic originating on Sq does not intersect Bq. We furthermore 
assume that {T>,g) may be isometrically embedded in some larger open manifold {V,g) 
such that g is a smooth extension of g across the boundaries ofV. Finally, we assume 
that for r > R the surfaces of constant r are timelike. See Figure [7| 
(ii) In the coordinates (r, t, the metric coefficients, assumed smooth, have the fol- 

lowing asymptotic behaviour (A,B £ {x^,y^}) 

gu = -i-l + 0{l), 5„ = g - 1^ + o (^) , 

(2-1) 9AB=r'ysAAB + 0{l), 

gt, = 0{i,), g,A = 0{l), g,^ = 0{^,). 

Here gfg2 is the standard round metric on the two sphere. It's convenient to also state 
the inverse metric components: 

gtt = _!^ + ^ + 0{^), <^- = r^ + i + 0(i), 
(2-2) g^^ = ^ Ws^f^ + 0{^), 

9''^ = 0{^), g'^ = 0{^), grA = o{^). 
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We finally assume that each radial derivative decreases the power of r by one. 

We make here a few remarks about these assumptions 

1. When the scalar field is coupled to the gravitational field assuming Neumann boundary 
conditions, the fall-off assumed above is too stringent. Fortunately, the fall-off assumed 
on the metric is actually far stronger than required for the theorems we shall prove. We 
could in fact assume 

9u = -'i + Oir), g,^ = J^ + 0{^), 
9AB = ^sAab + 0{r), 

5tr = 0(^), 9tA = 0{r), g,^ = 0(i). 

and our results will still hold. 

2. We do not need to make the topological assumption that X is cylindrical. We could 
equally well consider a conformal boundary with more elaborate topology or geometry. 

3. We have restricted to the case of four spacetime dimensions. Again, the theorems we 
shall subsequently prove do not require this restriction, so the results will apply in higher 
dimensions, with appropriately modified fall-off for the metric components and the field. 

4. If we have a spacetime for which the wave equation is well posed away from an asymp- 
totically AdS region (for example the AdS Schwarzschild spacetime), we can combine 
well posedness for the asymptotically AdS region with well posedness for the interior 
using the domain of dependence property of hyperbolic equations to show well posedness 
throughout. For more details, see the discussion in §8 of [9J. 

3. The modified energy 

We will now consider for a moment the case of the exact AdS spacetime. The usual 
energy one associates to solutions of the massive wave equation is given by 



(3.1) E[^t] = / T^uK'^dS'^, 

where is the timelike Killing vector and the energy-momentum tensor is given by 

(3.2) T^^ = V^0V^0 - ^g^u (Va(t>V^^ + ^0^^ , 



and satisfies VnT^'^ = when </> is a solution of ( 1 . 1 ) . If (/> has Dirichlet decay then one 



expects, by power counting, that -^[St] will be finite. However, if we have Neumann decay 



then the integral in (3.1 ) fails to converge near infinity. In order to deal with this problem. 



Breitenlohner and Freedman modified the energy momentum tensor (3.2) to give a new 
tensor 



(3.3) Tfj,,y = T^y + AT^,, = T^y + K{g^yU - V^Vy + R^i,) 

The modification satisfies 



k2 
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and so for the exact AdS spacetime, T^i^ will also give rise to a formally conserved energy. 
It transpires that the new energy differs from the original energy by a surface 

term. This surface term vanishes when (p decays like the Dirichlet branch, but diverges 
when (j) decays like the Neumann branch. By choosing k appropriately, it is possible to 
construct an energy which is positive, finite and conserved for both Dirichlet and Neumann 
decay conditions. 

Our proofs will make use of the modified energy associated to the timelike vector dt 
which, however, will only be approximately conserved since we no longer assume an ex- 



actly stationary spacetime. Rather than using the definition (3.3), it is in practice more 
straightforward to work directly from the PDE. In order to see how this works, we will 
briefly discuss a toy model which captures the salient features of the problem. 

3.1. A toy model. Consider the wave equation 

(3.4) utf-u^^- — + 0^4 = 0, 0<x<l, 

where < a < 1, subject to initial conditions 

n(2;,0) = no(x) 

(3.5) ut{x,^) = ui{x). 

Considering the behaviour near x = 0, we hope to impose as boundary conditions either 



u ^ X 



[l + 0{x^)] (Dirichlet) 
[l + 0{x'^)] (Neumann). 



At X = 1 we will require that n = 0. 

Suppose we have a suitably smooth solution to this equation. We can multiply (3.4) by 
xut and, after integrating by parts, deduce the conservation law for the standard energy: 

(3.6) ^ ^ Jt\ Jo (^^t + '>^l + a'^^^ ^dx = [xu^utll . 

For Dirichlet boundary conditions at x = the right hand side vanishes, however, for 
Neumann boundary conditions it is infinite. In order to introduce the modified energy, it 
is convenient to re-write the equation in the following form 

(3.7) u,,-x-i+"^(xi-2"^(x"n))=0. 



Which gives (3.4) upon expanding using Leibniz rule. We can multiply (3.7) by xut and 



integrate by parts to deduce 

(3.8) '^ = ±^.j\u] + [x-^d,{x^u)f) xdx = [xutx-°a,(x"u)]; . 







Now, for Dirichlet conditions at x = we again find that the right hand side vanishes, 
however we now find that it also vanishes for the Neumann behaviour at x = 0. The two 
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energies differ by a surface term: 

E- E = - \au 



If 211 

2 K Jo' 



wliich vanishes for Dirichlet conditions at x = and is infinite for the Neumann conditions. 
In this sense we can view as a 'renormahzed' energy, since we have formahy subtracted 
an infinite boundary term from the infinite energy to get a finite result. 

Thus even though the standard energy is infinite for the Neumann behaviour, we can 
modify it to get a conserved, finite, positive energy. We see now the justification for using 
the terms 'Dirichlet' and 'Neumann' to describe the boundary conditions. The Dirichlet 
condition requires u — )• as a; — t- 0, while the Neumann condition requires — t- 
as a; — )■ 0. 

This discussion also suggests that it will be fruitful to re-formulate the equation in terms 
of 'twisted' derivatives of the form x~"dx{x'^-). We shall do so in the next section and this 
will lead us to the appropriate setting in which to discuss the well posedness of (1.1). 



4. Well Posedness of the Weak Formulation 

4.1. Defining the problem. Motivated by the discussion of the previous section, we can 
now define the framework in which we shall work. 

We assume that U C is a bounded subset of with compact C°° boundary 
dU. This means that in the neighbourhood of any point P £ dU, there exists an open 
neighbourhood Wp C U of P and a smooth bijection $p : Wp — t- n B{0,6p), where 
= {(x, x") -.xy 0} and B (x, r) is the open ball centred at x with radius r. 

We're also going to assume that there exists a smooth function p : U ^ M+, which 
vanishes only on dU and such that there exists e so that if d{x, dU) < e, we have p{x) = 
d{x, dU) and if dU) > e, p{x) > e. We will set rii = dip, which extends the unit normal 
of dU into the interior of U . We may assume that the neighbourhoods Wp are such that 
po $p^(x,x") = x and 5p = e. We denote by Ut the timelike cylinder (0,T) x U and by 
dUr the boundary (0, T) x dU. 

We define our twisted derivatives in a similar vein to above. For a differentiable function, 
we set 

We may now define the equation in which we are interested: 
(4.1) 
wher^l] 

Cu = 

subject to the initial conditions 

(4.2) u{x,0) = uq, ut{x,0) = ui. 



utt + Cu = f in Ut, 
-dl'" {aijdfu) + bid"u + cu 



We use the Einstein summation convention, so that repeated indices i,j etc. should be summed over. 
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We assume all coefficients aij,bi,c are in C°°{Ut), however this is certainly stronger than 
necessary. We will assume throughout that aij is a symmetric matrix such that the uniform 
ellipticity condition holds: 

(4.3) e\C\^<a,,e^' 

for any ^'^ G R-^, where 9 is uniform in both time and space coordinates, and furthermore 
that aij is independent of t, on the boundary dU. 

For the time being, there are two possible boundary conditions in which we shall be 
interested. We will consider both Dirichlet: 

(4.4) n = on dUr, 
and Neumann: 

(4.5) nittijdj'u = on OUt, 

boundary conditions. 

To justify considering this equation, we have the following Lemma 

Lemma 4.1.1. Suppose iT^^g) is an asymptotically AdS spacetime patch with cosmological 
constant A = —p- Then there exists a smooth Lorentzian metric , g, on the solid cylinder 
Ut = [0, T] xi?(0, 1) CM.^ such that iV-,g) embeds isometrically into {UT,g), with I mapped 
to the boundary of the cylinder. Furthermore, setting <j) = -^^j^u for some p G C^{Ut) 
with p — >• 1 on the boundary, the wave equation 

(4.6) □^(A-^<^ = 



may be cast in the form (4-1) for some p, Oij, hi, c satisfying the assumptions above, with 



p = r 



-^ + 0(r-^) anda = a/| + A. 



Note that for — | < A < — | we have < a < 1. Making use of this lemma, we see that 



a well posedness result for (4.1) implies a well posedness result for the asymptotically AdS 
space by restricting to data supported sufficiently close to the boundary, where we assume 
trivial data on Bq. 

4.2. The function spaces. In order to introduce a weak formulation for the initial- 
boundary value problem we are considering, it will be necessary to define the function 
spaces in which we seek a solution. For a locally measurable function u and measurable 
set V C U, we define the norm and space 

(4.7) \\u\\l2^y^ = j^u'^pdv. L^(y) = {u : ||u||^2(y) < oo}. 
This is clearly a Hilbert space with inner product 

(4.8) (ui,U2)iJHy)= j uiU2pdv. 
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Now, we note that for smooth functions (j), ip of compact support we may integrate by parts 
to find 

/ (t)dfi}pdx = - [ {d]-'^<t))i^pdx, i = l,...,N 
Jv Jv 

This allows us to define a weak version of d". We say that Vi = d^u is the weak a- twisted 
derivative of u if 



(4.9) / Vi(j)pdx = - I udl~°'(ppdx 



for all (j) G C^{V). We say that u G H}{V) if d°u exists in a weak sense and dfu G L^{V). 
We can define a norm and inner product on H^{V) as follows: 



N 
i=l 

N 



(4.10) {uI,U2)h\V) = (^^l'^^2)L2(y) +X](^f^l'^i''"2)L2(y) • 



Next we define H^{V) to be the completion of C'^{V) with respect to the norm |M|//i(y)- 
We shall often take V = U . On any subset compactly contained in [/, these spaces are 
simply equivalent to the standard Sobolev spaces. 

We note at this stage that dl~'^ is the formal adjoint of df with respect to the inner 

product. Thus the second order operator dl~" (^aijdj-^ appearing in (4.1) is formally 
self-adjoint. When we come to consider higher regularity, we shall need the Sobolev space 
associated to the adjoint derivative operator. In particular u G H_ (V) if dl~'^u exist in 
a weak sense and dl~°'u G Ll^Y)-, with the obvious inner product and norm. Again we 
define ^o(^) be the completion of C'^{V) with respect to the norm 
Let us state some properties of functions in these spaces. 

Lemma 4.2.1. (i) Functions of the form u = p~"v, with v G C°°{U) are dense in 
H\U). 

(a) Suppose u G H}j{U), then d^u G Ll{U) for any j3. 

(Hi) Suppose u G H}{U). Then in a collar neighbourhood [0,e) x dU C U of the boundary, 
we have u G C*^((0, e); L^(5C/)), with the expansion 

u = p~''[uo + 0(p'')) 

where uq G L?'{dU), with uo = if u € H}^{U). Furthermore, for any 5 > 0, there 
exists a Cs such that 

(4.11) WuoWmdU) < ^M\h^u) + <^s\H\lHu) 

Similar results hold for H_, but with a replaced by 1 — a. 
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Part (i) follows from a result of Kufner [13j, and parts may be derived by 

showing that the inequalities hold on suitable dense subsets. From this we see that if 
u G H}{U), then u may 'blow up like /9~" near dU\ whereas if u S H}){U) then u is 
'bounded near dU^ in some appropriate sense. These spaces thus capture, to a certain 
degree, the boundary behaviour we hope for in our solutions. A consequence of the proof 

of (ii) is that ^o(^) = kI{U). 

In fact, we can prove a sharper result about the range of the trace operator, together 
with an extension result: 

Lemma 4.2.2. The operator T o , where T is the trace operator, maps H^{U) into 
H"(dU), and the map is surjective. Furthermore there exists a bounded right inverse so 
that corresponding to any uq G H°'{dU), there exists a u £ H}{U) with p'^u\qu = uq in the 
trace sense, with the estimate 

\\'^\\h^{u) - ^W'^oWn^'idU) 

where C is independent of uq. 

This follows from the fact that p^'u belongs to a weighted Sobolev space, to which one 
may apply the results of [Tl] . 

We will also require the spaces H^{U)* and H}{U)*, the dual spaces of H_q{U) and 
H^{U) respectively. If f eX*,ueX we denote the pairing by 

{f,n), 

and define 

11/11^. = sup{(/,u) \u G X, \\u\\x < 1}. 

It will be convenient, for notational compactness, to define the following spaces and 
norms. The Neumann data space -f^^^ta n(^) consists of triples (uq, ui, f) with uq G H^{V), 
ui G L^{V) and f G L^{[0,T]; L^{V)), whereas for the Dirichlet data space -ffjata n(^) 
additionally require uq G H^q. For both spaces, we define 

||(uo,Ul,f)||^l^^^(^) = ||U0|||-l(y) + ||ltl|||2(y) + l|f|li2([o,T];L2(y)) , 

We take H^^^^ ^{V) to consist of u G L'^{[0,T]; Hl{V)) with 

Illl2 ||||2 ||||2 ||||2 

"""^iol., d(^) ^ ll"llL°°([0,T];/i;i(y)) + II"IIw^i>°°([0,T];L2(F)) + W^W {[0,T];{H},(V))*) < °°' 

and ffgQi j^{V) to consist of u with 



'""^sol. n(^) ~ ll"lli°°([0,r];Hi{y)) + ll"lll4^i.°°([0,T];L2(y)) + \\^\\H2{lO,T];(H}iV))*) < 



4.3. The Weak Formulations. In order to motivate the definition of the weak solutions, 
let us suppose that we have a solution to 

(4.12) uu + Cu = f onU 
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which is sufficiently smooth for the following operations to make sense. We can multiply 
the equation by a smooth function f , integrate over U and integrate by parts to establish 

Ju Ju JdU 

The surface term will vanish either if u satisfies the Neumann boundary conditions, or else 
if V satisfies the Dirichlet conditions. We define the following bilinear form on H^(V) 

(4.13) Bv[u,v;t] = [ [aij{d^u)idfv) + b\d^u)v + cuv] p dx. 

Jv 

If B has no subscript, we assume the range to be U. Now we may define the weak Dirichlet 
and Neumann problems: 

Definition 2 (Weak Dirichlet IBVP). Suppose {uQ,ui,f) G H^ata dO^)- ^"-V ^^^^ 

u G Hlgi j:){U) is a weak solution of the Dirichlet IBVP if 

i) For all v G H_o{U) and a.e. time < t < T we have 

(ii, v) + B[u, v; t] = (f , v)^2f^^^ . 

a) We have the initial conditions 

u(0) = 'Uo, u(0) = ui. 

Definition 3 (Weak Neumann IBVP). Suppose (uo,ui,f) G H^^^^ ^i^)- ^^^^ 
u G H'^^i jv(^) ^ weak solution of the Neumann IBVP if 

i) For all v G Il}{U) and a.e. time < t < T we have 

(ii, v) + B[u, v; t] = (f , v)^2fjj^ . 

ii) We have the initial conditions 

u(0) = uo, u(0) = ui. 

We note that by the calculation above, a strong solution obeying the Dirichlet (resp. 
Neumann) condition on the boundary is necessarily a weak Dirichlet (resp. Neumann) 
solution. The converse of course need not be true, however if we have enoug h reg ularity 
to integrate by parts then taking an arbitrary v G II_o{U) we conclude that (4.12) holds 
almost everywhere in U in both the Dirichlet and Neumann case. Noting in the Neumann 
case that the trace of is arbitrary on the boundary we can, with care, deduce that 
UiOijd^u G Hq{U). We will see this in more detail later when we consider the asymptotics 
of the solutions. 



4.4. The Weak Formulations. 
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4.5. The theorems. We're now ready to prove the well posedness of solutions to the weak 
formulations of (4.1). First, we have the following result 

Theorem 4.1 (Uniqueness of weak solutions). Suppose u is a weak solution of either the 
Dirichlet IB VP or of the Neumann IB VP. Then u is unique. 

Proof. The proof of uniqueness for the weak solutions procedes almost identically to the 
proof of uniqueness of weak solutions to a finite IBVP. Without loss of generality, one may 
assume trivial data. In both cases one may take as test function 

("^) v,).{/>M- »|;|». 

and then integrate the weak equation over < t < s. Standard manipulations making use 
of the uniform hyperbolicity condition (4.3) then show u = 0. For example, one may take 
the proof of Evans [111 p. 385] and replace the standard spatial derivatives with twisted 
derivatives. □ 

Theorem 4.2 (Existence of weak solutions). (i) Given (no,ni,f) G Il\ata d(^)' ^^^'^^ 
exists a weak solution to the Dirichlet IBVP corresponding to this data, 
(a) Given data (no,^ii,f) G -f^^ata Af(^)' ihere exists a weak solution to the Neumann 
IBVP corresponding to this data. 
In both cases, we have the following estimate 



where C depends on T, [/, a and the coefficients of the equation. ] stands for D or N as 
appropriate. 



It is convenient to divide the proof of Theorem 4.2 into several Lemmas. We start by 
picking a sequence < < e which decreases monotonically to zero, and define the sets 
Vk = {x : p{x) > Ofc}. The broad strategy is to solve the finite IBVP on each 14, where 
the equation becomes strictly hyperbolic and classical theory applies, and then find a way 
of passing to the limit 'A; — )• oo'. 

Lemma 4.5.1. (i) Data (uq, ui, f) such that for all k > k^ the problem 

utt + Cu = f on [0, T] X Vfc, 

(4.15) u = uo\y^ ,ut = uo\y^ u = on {0} x 14, u = on [0,T] x dVk 

has a solution which is C°°([0,T] x Vk) form a dense linear subspace of H^^^^ _d(^) 
(a) Data (uoi ^^ii /) such that for all k > ko the problem 

Utt + Cu = f on [0, T] X Vk, 

(4.16) u = uo\y^ ,ut = uo\vk u = on {0} x Vk, UiOijd'^u = on [0,T] x dVk 

has a solution which is C°°([0,r] x V^) form a dense linear .subspace of H\^^^ n(^) 
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Proof. For (i), we may take uo,ui £ C^{U) and f G C^{[0,T] x U). For large enough 
k the data are supported inside Vk, and data of this form a dense Unear subspace of 
-^data d(^)- (^^)' ^® need the fact that smooth functions u for which riiaijdj'u = 
outside a compact set are dense in H^{U). We take uq to have this property, and take 
ui G C^{U). FinaUy we can take a smooth f which is a sum of one component in C^{U) 
and another of arbitrarily small L^([0, T]; L'^({7)) norm which ensures the higher order 
compatibility conditions vanish to all orders on t = 0. For large enough k this data will 
launch a smooth solution and such data are dense in H^^^^ n(^)- '-' 

Lemma 4.5.2. Suppose u is a smooth solution of 

u'[t + Cu'' = f on[0,T]xVk, 
^ - uo\y^ , = ui\y^ on {0} X Vk, 



w 



(4.17) 

with either = or UiOijd'ju'^ = on dVk- Then satisfies the estimate 



(4.18) 



u 



<C||(uo,ni,f)||^i 



where C is uniform in k. \ stands for D or N as appropriate. 

Proof. Multiplying by and integrating by parts, using the boundary condition to neglect 
the boundary term, one has 



(4.19) 

We also note that 
(4.20) 



+sv'Ju^u^^] = (f,u 



d_ 
dt 



u 



2 



I?{U) 



We add these two equations and make standard manipulations, using the uniform hyper- 
bolicity condition, together with Gronwall's lemma as in the case of the finite IBVP (see 
im p 381.]) to derive the estimate (4.18). Careful examination of the steps shows that 
one may take C uniform in A;. □ 

Lemma 4.5.3 (Weak compactness). (i) Suppose G ^lo\T)i^k) , with 

< C 



u 



Then there exists u G -f^soiD(^) ^^^^ 1 1^1 Ih^ j ^(i/) — ^ ^"-^ ^ subsequence u^'- such 
that for any v G ^o(^), taking I large enough that ki > m we have for almost every 
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t: 

(4.21) (u'''{t),v 

u''^{t),v 

(a) Suppose u'' G -f^sol.N(^fc)' "^^^^ 



u 



(u(t),z;)^2(^) 
(9fu(t),t;)^2(j;) 

(U(t),f)^2(^) 

(u(t),?;), 
< C 



Then there exists u G -f^soiN(^) ^^^^ ll^llf/i ('(ti — ^ ^'^'^ ^ subsequence u^' suc/i 

sol.N^ '' 

that for any v G H}{U), we have for almost every t: 

{dfu{t),v)jj,^u) 

(U(t),u)^2(j;) 

{uit),v), 



dru'^it),v\y^ 



(4.22) 



Proof We demonstrate first the proof for G :H^^(Vfc), ||ii||iifi(yj.) < C, i.e. the first part 
of (ii). We define ^ €L^{U) to agree with on Vk and to vanish on U\Vk. Similarly, we 
define d^u'' G L^{U) to agree with dfu'^ on Vk and to vanish on U\Vk. Weak compactness 
of L^{U) gives a weakly convergent subsequence [u^^ ,dfu^i) {u,Vi). It remains to show 
that Vi = dfu in the weak sense. To show this, multiply dfu^^ hy (j) & C^{U) and integrate 
over U. For / large enough that supp C Vfc, , we have (pd^u^^pdx = — dl~'^4>u^' pdx, 
so by taking weak limits we're done. Similar considerations may be applied to the other 
results in the Lemma, after applying Riesz representation theorem to ii. □ 

Remark : This Lemma can be extended to apply to higher spatial derivatives of in 
an essentially unchanged fashion. 

Now we can combine the results above to show that there exists a solution to the weak 
problems. 



Proof of Theorem ^.2. (i) Suppose we have data {uo,ui,f) such that for all k > ko the 
problem 

utt + Cu = f on [0,r] X Vk, 
(4.23) u = uo\y^ ,ut = uqIv,. w = on {0} x V^, u = on [0,T] x dVk 



THE MASSIVE WAVE EQUATION IN ASYMPTOTICALLY ADS SPACETIMES 

has a solution, u'^ which is C°°([0,T] x Vk)- By Lemma 



4.5.2 
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we have the estimate 



u 



< C||(U0,Ml,f)||j:^l^^J^) 

And we also know that for k > m and for any v G H_Q{V.m) we have 
(4.24) (n,^„ v) + By, [n^ v; t] = (/, v)^2^y^) 



Applying Lemma 4.5.3 we conclude the existence of u satisfying 



(4.25) 



u 



, ^(m < C||(uo,ui,f) 

and for any v G H_Q{Vm): 
(4.26) {utt,v) + B[u, V- 1] = (/, v)^2^jj^ . 

Noting that functions v E H^{Vm) are dense in H^{U), we conclude that u satisfies 
the first condition to be a weak solution of the Dirichlet IBVP. We must now check 
that the weak solution we have constructed satisfies the initial conditions. Fo r this , 
choose any function v G (7^(0, T; C^(f/)), with v(r) = v(r) = 0. Integrating ( |4.26 ) 
in time, we have after twice integrating by parts 



(i;,u) + B[a,w]t\dt 



dt-(u(0),v(0)) + (u(0),v(0)) 



similarly, we have from (4.24) 



i (^'U*^)^,^+^vJu^v;^]d^ = ^ (f,v)^2(^,^)dt 
Setting k = ki, passing to the limit we have: 



u'^(0),v(0)) +(u'^(0),v(0) 



^ {v,u) + B[u,v;t]dt = (f,v)^2(^)dt-(M0,v(0)) + (txi,v(0)). 

Since v(0), v(0) are arbitrary, we conclude that u(0) = uq, u(0) = ui and we're done. 
Finally, we make use of Lemma 4.5.1 together with the uniqueness result Theorem 



4.1 and a standard argument based on continuity, using (4.25 ), to show that our result 



holds for any (no,ni,f) G H^^^^^ ^{U). 
(ii) The Neumann case follows in an almost identical manner, solving a sequence of finite 
Neumann problems for suitably smooth data and using the weak compactness to 
extract a weak solution. Finally a continuity argument again extends the existence 
proof to all admissible data. 

□ 

5. Higher Regularity and Asymptotics 

We now wish to show that if more assumptions are made on the data, the weak solution 
can be shown to have improved regularity. In order to do this, we will require some elliptic 
estimates, enabling us to control some appropriate norm of u in terms of Cu. 
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5.1. The norm. The norm is sHghtly unusual in its definition because it is neces- 
sary to distinguish the directions tangent to the boundary from those normal to it. We fix 
a finite set of vector fields {T^^^, A^^^^j on U which satisfy the following properties: 

(i) For p < €, we have T^^^ normal to dip, while N^^^ are parallej^to Uijdjp. 

(ii) At each point of 17, the set {T^^\n'^^^ spans M^. 

Definition 4. We say that a function u G H}{U) belongs to H^{U), provided 
for all A, B, and we define the norm: 



h\u) 



Remarks : 

(a) A different choice of {T^^^^, A^^^^} satisfying (i), (ii) gives rise to an equivalent norm. 

(b) If n G H^{U) then u G Hl^ {U). 

(c) If n G H^{U), then G iJ^iU)- This observation is important in 
establishing the higher regularity energy estimates. 

5.2. Elliptic estimates. We first define the weak version of the elliptical problem we 
study. We assume that t is fixed throughout this section: 

Definition 5. Suppose f G H^{U)* (resp. H}{U)*). We say that u G H}i{U) (resp. 
H}{U)) is a weak solution of the Dirichlet (resp. Neumann) problem 

(5.1) Cu = f onU 
with u = (resp. UiOijdj'u = 0) on dU, if 

B[u,v] = {f,v) 

for all V G HI{U) (resp. H^{U)). 

Theorem 5.1 (Ellipt ic E stimates). Suppose u is a weak solution of either the Dirichlet 
or Neumann problem (5.1) and suppose that in fact f G L^{U). Then u G H(^{U) with the 
estimate 

(5.2) \\u\\HHm<c{\\f\\LHu) + \W 



Furthermore, in the Dirichlet case Tidfu G H}^{U) and in the Neumann case Nidf G 

t;(f/). 

We split the result into several Lemmas 



Lemma 5.2.1. There exist constants Ci,C2 and //o > such that 
(i) \B[u,v]\ < Ci ||n|||^i(^) 

^recall dij = on dU, so this condition is independent of t 
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(a) C2 ll-ulll-i^f;) < B[u,u] +/io||^^|||2(^) 

Proof. This is a standard manipulation, making use of the uniform ellipticity of ajj. □ 
Lemma 5.2.2. There exists /xq G M such that for all fi > fio, f G l/{U) the equation 

Cu + iJ,u = f, on U 

with either Dirichlet or Neumann boundary conditions, has a unique weak solution satis- 
fying 

\\u\\h\u) <C\\f\\L2^u) ■ 

Proof. Because of the estimates in the previous lemma, we may apply the Lax-Milgram 
theorem to Bu[u, v]+ {u, v)^2fjj^ thought of as a bilinear form on either H^{U) or H}{U) 
for Dirichlet, Neumann conditions respectively. □ 

Lemma 5.2.3. Suppose f £ C^{U), and k is sufficiently large that supp f C Vk- Then 
for n> Ho 

(5.3) Cu + = f, on Vk 

has a unique solution in C°°{Vk). Furthermore, this solution obeys the estimate 

(5-4) MHHV,)<C\\f\\L2^u) 
with C uniform in k. 

Proof. We can apply Lemmas 5.2.1, 5.2.2| on Vk to deduce the existence of a unique weak 
solution to (5.3) with the appropriate boundary conditions, satisfying 



\\u\\h\Vu) <C'II/IIl2(c/)- 

where C is independent of k. On Vk, the operator C is uniformly elliptic in the standard 
sense, so classical elliptic estimates imply that u is smooth. We have that 



where 



2 < C' ||/||j;^2(^^. Focusing on a coordinate patch, we can work assuming 

Ll (^fe) ^ 

= {{x,x"') : X > ek, + x"'x'^ < e}, p = x and assume that C is a smooth cut-off 
function on {(x,x") : x > 0, + x'^x'^ < e} which vanishes on the curved part of the 
boundary. We note that 9" = da and that d"da = dadf. Now consider the following 
integral, where the index A S {2, 3, . . . , A^} is a fixed index, with no summation over it. 

(5.5) I ■= I [dl-''{aiMu)][dA{C'dAu)]xdx 

< 5\\C'^dAdAu\\^2^y^ +C (\\f\\l2^v) + \Mh}{v)^ 

By choosing C large enough, we may take 5 to be arbitrarily small. Now, since C vanishes 
on the curved part of dV and either u or niOijd'^u vanishes on the flat part, we may 
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integrate by parts twice to find 



(5.6) 



> 



[dA{aijdfu)][d2{edAu)]xdx 



> e\\C^d^dAU 



C\\u\ 



where in the last hne we have used the uniform ehipticity of Uij. The constant C here 
depends on the functions aij,C,, which are uniformly bounded in k. Now taking (5.4), 



( |5.5[ ), (5.6) together, and choosing 6 sufficiently small, we have that 
(5-7) ||C'5a^||^1(v^^)<C||/||^.(^) 

with C uniform in k, so we have estimated the tangential derivatives. Returning now to 
the equation, we can write 

^l-''{a.^^tu) =f- ^a{aa^^fu) 

Multiplying by C^, we readily estimate 

(5-8) ||C'a„9>||^i^^^^<C||/||^2(^). 

Combining these estimates with a partition of unity subordinate to a set of coordinate 
patches covering the boundary and an interior estimate which follows from standard elliptic 
theory, we're done. □ 



Proof of Theorem 5. 1 . First we note that if u is a weak solution of ( 5. 1 ) with either Dirichlet 
or Neumann boundary conditions, then u is the unique weak solution of 

(5.9) Cu + fiu = f 

with / = / + ^M for sufficiently large /i. Suppose / G C^{U). Then we can solve the finite 
problems on Vk, with the estimate 



< C 



with C uniform in k. We deduce the existe nce of a subsequence which tends weakly to u 
in H^, as in the remark after Lemma 4.5.3, so we find that 



< c 



f 



I?{U) 



now, we may relax the condition that / G C^{U) since such functions are dense in L^{U). 
Finally, replacing /, we deduce 



< C 



Finally, note that for Dirichlet conditions we have G Hs^iVk) and for Neumann we have 
niaijdj'u^ G H_l{Vk), so that in the limit u G H_q{U) or niUijdj'u G H^{U) respectively. □ 
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We would like to also prove elliptic estimates at a higher level of regularity than H^. 
Unfortunately, the behaviour of the solutions near the boundary doesn't lend itself to a 
description in terms of a global Sobolev space. Accordingly then, we first consider interior 
regularity. 

The orem 5.2. Suppose u is a weak solution of either the Dirichlet or Neumann problem 
{5l\ l and suppose that in fact f e L^{U) n H^^^XU). Then u G H^{U) n H'^^'^{U) 

Proof. This follows from standard elliptic estimates and the fact that C is uniformly elliptic 
on any V CC U . □ 

To say more about the behaviour near the boundary, we shall once again need to dis- 
tinguish directions tangent and normal to the boundary. It's convenient to introduce the 
space H]p{U), consisting of all functions u such that 

r«r(2)...r«nGL2(c/) 

for any / < m smooth vector fields T^*-* tangent to the boundary dU. To capture the 
behaviour normal to the boundary, we work with an asymptotic expansion. 

Theorem 5.3. Suppose u is a weak solution of either the Dirichlet or Neumann problem 
(5.1) where f G H^(U) D (U), where m > 0. Suppose further that if m > 1 near dU 
we have the following expansion for f: 

(5.10) / = + + + V+_i + 0(p'"-"))] 

[fo + P/r + • • • + P'"~'/^-i + O (p"-^+"))] 

where 

ff £ H'^-^-'idu) 

and\\p-''0{p'')\\j^2(^QU) is bounded as p^O. Then u e (U) n R^^^^ (U) n H^;^^ {U) has 
the following expansion for m > 0: 

(5.11) u = p"[n+ + pn+ + ... + p-n++i + 0(p™+^-))] 

+p~^ K" + pnr + . . . + p'^'-'u-^, + O )] 

where 

uf £ H'^+^-'idU). 

Furthermore if u satisfies the Dirichlet conditions, Uq = = 0, while if u satisfies the 
Neumann conditions = 0. 



Proof. The proof is by induction. To establish the m = case, we apply Lemma 4.2.1 and 
Theorem 5.1 to deduce that in a coordinate patch near the boundary 

a^idfu = 2;-"[0(x")]+2;"-i[c+ + 0(xi~")], 

(5.12) daU = j;-°[c_+0(x")]+x"-^[0(x^-")]. 

with c± G Lp'{dU). We thus have that 
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and integrating this gives (5.11) for m = 0, with uf,UQ,uY G L'^{dU). Finally we note 
that the second identity of (5.12) implies Mq G H^{dU). In order to get the induction 
step, we first commute with a vector field tangent to the boundary, which establishes all 
but the highest order in p of (5.11 ) by the induction assumption. To get the highest order 
terms, we re-arrange the equation Cu = / to give an equation for 5j.~"(9^M, making use 
of the induction assumptions and integrate twice. Taking care of the boundary conditions 

= 0, while for Neumann 

□ 



imposed shows that for Dirichlet conditions, we have Uq 



0. 



Taking a little more care about the origin of terms in the series, we can easily show 



Corollary 5.4. (i) Ifu, f satisfy the conditions for Theorem 5.3 with Dirichlet boundary 
conditions and furthermore f^ = for < z < m — 1, then = for < i < m + 1 . 
(a) Ifu, f satisfy the conditions for Theorem 5.3 with Neumann boundary conditions and 
furthermore f^ = for < i < m — 1, then uf=0 for 1 < i < m + 1. 

5.3. Higher regularity. We define the higher regularity data spaces inductively as fol- 
lows. We say (no,ni,f) G Hj^^^ j^ if S HJ{U),ui £ Hl{U),i £ H\[0,T]; l/{U)) with 
the product norm. In the Neumann case, (uo,ui,f) G -ffJata at if ^o S H^{U),aijd'^UQ G 
Hl{U),ui G W-{U),i£ H\[0,T];L^{U)). Next, we define 

50 = Uo gi= ui 



(5.13) 



9i+2 



E 

1=1 



Here is the second order operator given by differentiating the coefficients of C i times 
with respect to t. For m > 2, we say (uo,ni,f) G H^a.ta., f (^o,wi,f) G H^^'^ p fW G 
L2([0,r];i/™7^-i(?7)) for < ^ < m - 1 and (5m-i, 5™, f^™"')) G ^La, f Here as usual 
t stands for D or N as appropriate. We define the norms 

2 



{uo,ui,i)\\fj 



data, t 



{uo,ui,i)\\^m-i 



+ 



data, t 



data, t 



these spaces are chosen so that the relevant 'compatibility conditions' hold. We may show 
that if {uo,ui,i) G F^^ta, t' then uq G KJU), m G H^'^iU). 

Theorem 5.5. (i) Suppose u is a weak solution of the Dirichlet IBVP corresponding to 
data (no,ni,f) G H]iata oi^)- Suppose in addition, {uo,ui,i) G H'^ata d(^) ^^^^ 

u G L^{Q,T-H\U)), u G L-(0, T; ^^(C/)), 
ii G L-(0,r;L2 ([/)), u G L^{0,T; (HKU))*), 

with the estimate 



(5.14) 



ess sup 

0<t<T 



|u(t)||j 



+ l|u(t)|| + l|u(t)|| 



+ u 



L^{0,T;iHliU)r) ^ C'll(^0,Wl,f)|lH2 



data, dC^) 
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(ii) Suppose u is a weak solution of the Neumann IBVP corresponding to data {uq, mi, f) G 
^\ata n(^)- Suppose in addition, (uo,ni,f) G H'j^^^ Ar(^) ^^^^^ 

u e L'^{0,T;H^{U)), u e L'^{0,T;H\U)), 
lie L-(0,T;L2(C/)), ueL\0,T;iH\U))*), 



with the estimate 



(5.15) 



ess sup 

0<t<T 



H'^(U) 



+ u 



+ u 



L2 



(0,T;(ifi({/))*) < C||(uo,Ui,f)||^2 ^(^) 



Furthermore UiOijOfu G L°°([0, T]; iio(C/)). 

Proof. First note that without loss of generaUty, we may take uq = 0, so that data which 
give rise to a smooth solution to the restricted problem on Vk for k sufficiently large 
are again dense. We return to the approximating sequence u'^ we established in proving 
the existence of a weak solution. Commuting the equation with dt and making use of 
the elliptical estimates on Vk established in the previous section it is straightforward to 



derive bounds for li 



lL°°{[0,T];|/i{yfc))' 



u 



\L^{[0,T];I?iVk)) 



and the relevant norm of u 



which are uniform in k. Passing to a weak limit and applying Theorem 5.1 to deduce 
u G L°^([0,r];F2(;7)), we're done. □ 

Commuting further with dt, it can be shown that the following theorem holds: 

Theorem 5.6 (Higher Regularity). (i) Assume (uo,Mi,f) G ^dataO ^^'^ suppose also u 
is the weak solution of the Dirichlet IBVP problem with this data. Then in fact 



(5.16) 
(5.17) 



ess sup 

0<t<T 




+ 



ir\u) 



m— 1 



U 



+ 



dt^ 



L\U), 



LHO,T;iHhm*) 



<C||(no,ni,f)||^^ 



data,D 



where C is a constant which depends on T and a and the coefficients of the equation. 
Furthermore 

u«(t) G ^rrr*([/) for 0<i<k. 



(ii) Assume {uo,ui,{) G HJ^ta N '^^'^ suppose also u is the weak solution of the Neumann 
IBVP problem with this data. Then in fact 



(5.18) 
(5.19) 



ess sup 

0<t<T 




+ 



d 



m—l 



U 



dt 



m—l 



+ 



H}{U) 



(i™u 



df"- 



L-^(Q,T;(H^{U)Y) 



< C\\{uo,Ul,i)\\fjr, 
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where C is a constant which depends on T and a and the coefficients of the equation. 
Furthermore 

u»(t) e H'^-'iU) for < i < m. 

Note that we do not directly control higher spatial derivatives of u in global Sobolev 
norms, although we do have control of powers of the C acting on u. We can however make 
use of (a very slight adaptation of) Theorem 5.3 to give the following asymptotic expansion 
in the situation where f has the appropriate behaviour near the boundary: 



Theorem 5.7. Suppose u satisfies the conditions of Theorem 5.6 (i) or (ii). Suppose that 
G L'^{[0,T];H^-'-^{U)nH'^-'-^{U)). Suppose further that if m > 3 near dU we have 
the following expansion for f : 

(5.20) / = + + + 3/+_3 + o(p— 2-))] 

[fo + p/r + • • • + p"-'/-_3 + o (p'"-^+"))] 

where 

Then u has the following expansion for m > 1.' 

(5.21) u = p"[^+ + pn+ + ... + p— \+„i + 0(p— ^-))] 

+p-- [u, + pnr + . . . + + O (p'"-i+"))] 

where 

Furthermore if u satisfies the Dirichlet conditions, Uq = = 0, while if u satisfies the 
Neumann conditions uf = 0. 

For a formal power series approach to determining the coefficients of these expansions, 
see ITS). 



6. Other boundary conditions 



We have now established well posedness and a regularity theory for solutions of (4.1) 
subject to either Dirichlet or Neumann homogeneous boundary conditions. We will discuss 
briefly some of the other possibilities listed in the introduction, although we shall not go 
into quite so much detail. 

6.1. Inhomogeneous boundary data. First, we define the weak formulations for the 
inhomogeneous problems. We assume throughout this section that uq G H}{U),ui E 
L2(C/),f G L2([0,r];L2(J7)). We furthermore take go,gi G i[0,T]; L^{dU)) to be some 
functions on the boundary. The gi will need to be subject to further conditions in order 
to give a well posed problem. 
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Definition 6 (Weak Inhomogeneous Dirichlet IB VP). Suppose uo,ui,{ are as above with 

the additional condition p"uo\gy^ — go(0). We say that u G L°°{[0,T]; H^{U)) with 
u G L°°{[0,T];L^{U)), ii € L'^{[0,T]; {HI{U))*) is a weak solution of the inhomogeneous 
Dirichlet IBVP: 

utt + JCu = f in Ut 

P'^'^\dUT = 30 on dUr 

u = uq, Ut = ui on {0} X U 

provided 

i) For all v G H^{U) and a.e. time < t < T we have 

(ii, v) + B[u, v; t] = (f, v)^2(^jj^ . 

ii) We have the initial conditions 

u(0) = no, u(0)='ui. 
Hi) We have the boundary condition 

p"^\dUT = go 

Note that in contrast to the homogeneous Dirichlet problem, the unrenormalized energy 
will be infinite for a solution of the inhomogeneous Dirichlet problem with go ^ 0. 

Definition 7 (Weak Inhomogeneous Neumann IBVP). Suppose uq, ui, f are as above. We 
say that u G L°^{[0,T]; H^{U)) with u G L'^{[0,T]; L^{U)), ii G L^{[0,T]; {H^{U)y) is a 
weak solution of the inhomogeneous Neumann IB VP: 

Utt + = f in Ut 

p^'~°'niaijd'^u = gi on dUx 

■' dUr 

u = no, Ut = ui on {0} x U 

provided 

i) For all v G H_^{U) and a.e. time < t < T we have 

{u,v) + B[u,v;t] = (f,u)i2(f;) + (gi, p'^v\g^)^2^Qu) . 

ii) We have the initial conditions 

u(0) = no, u(0) = ni. 

If we assume sufficient regularity, it is possible to show that the weak solutions are 
equivalent to strong solutions, by a standard integration by parts. 

It is clear that if ?; G L'^{[0,T]; H^{U))nH^{[0,T]- H^{U))nH^{[0,T]- I?{U)) = H^(Ut) 
satisfies the condition 

P°^v\qu = go for inhomogeneous Dirichlet, or 



p^-'^aijdfv 



= gi for inhomogeneous Neumann. 

dUr 
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then we can apply our previous weak well posedness results to the functions 

u = u — V 

We'd like to know what conditions are required on go,gi in order that such a v exists. The 
following Lemma gives the results we require, and comes from adapting Lemma 4.2.2 to 
H^{U). 



Lemma 6.1.1. (i) Suppose v G H^{Ut), then p'^v\qu and p 



trace sense, and we have 



dUr 



exist in a 



m+°'{dUT) 



+ 



3 \dUj 



m-<^{dUT) 



<C\\v\ 



KHUt) 



(a) Suppose vq G H^^°'{dUT) and vi G "^{dUx)- Then there exists v G H^{Ut) such 
that 



p'-'^aijdfvl 



Vl 



where the restriction is understood in the trace sense. Furthermore we may choose v 
such that 

\\v\\h^{Ut) - ^ (iboll/^i+c^Ot/r) + W'^'i-W m-^{dUT) 
with C independent ofvi. 

Taking this with our previous results, we conclude 

Theorem 6.1. (i) Given uo,ui,i as above, qq G H^^^'^dUx) with p'^uq\qu^ = 5o|f=ov 
there exists a unique weak solution to the Dirichlet IB VP corresponding to this data, 
with the estimate 



(a) Given uo,ui,i as above, gi G H'^~°'{dUT), there exists a unique weak solution to the 
Neumann IB VP corresponding to this data with the estimate 



We shall not go through the proof in detail, but it is clear that the higher regularity 
results of ^ can be extended to the inhomogeneous case. 

6.2. Robin boundary condition. Another possibility for a well posed boundary condi- 
tion is what we might call a Robin boundary condition: 

p^-"n,a„-9> + Pp^u = OondU 



(6.1) 



where f3 is some suitable function on dUx- We assume (3 G C^^dUx) to be concrete, but 
this is not necessary. This can be achieved in the weak formulation in a similar fashion to 
the introduction of an inhomogeneity for the Neumann condition. 
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Definition 8 (Weak Inhomogeneous Robin IBVP). Suppose uo,ui,i are as in ^6.1\ We 

say that u G L°°([0, T]; ^^(C/)) with u G L'=^{[0,T]; L^{U)), u G L^{[0,T]; {H}{U))*) is a 
weak solution of the Robin IBVP: 

Utt + = f in Ut 

p^^'^niaijdfu + Bp'^u] =0 on OUt 

^ J dUr 

u = Moi Ut = ui on {0} X U 

provided 

i) For all v G H}{U) and a.e. time < t < T we have 

{u,v) + B[u,v;t] + {p°'u\q^ , p°'v\g^)^2(^Qu) = i^,v)^2^jj) . 

ii) We have the initial conditions 

u(0) = uq, u(0) = Ul. 

It is straightforward to show that the well posedness and regularity results of ^ and ^ 
can be extended, where we require the estimate (4.11 ) to deal with the surface terms which 
arise. 

Remark: In the case where 6j = 0, c > and Ojj,c independent of time, we can relate 
our result to the theory of essentially self-adjoint operators. A consequence of Stone's 
Theorem (see for example [16]) for self-adjoint operators states: 

Theorem 6.2. Let L : D ^ H be a densely defined positive symmetric operator. Suppose 
for every f,g & D there exists a twice continuously differentiable solution (in D) to 

Utt + Lu = 0; u{0) = f; ut{0) = g 

Then L is essentially self-adjoint. 

In our case H = L^{U), L = C, however this operator is not essentially self-adjoint on 
C^{U). Thus, the choice of dense subspace D determines a self-adjoint extension of C. If 
we take 

D = {u = p^'v: V e C^{U),p'^d'u = OondU for k = l,2,.. .}, 

this gives the self-adjoint extension corresponding to Dirichlet boundary conditions, whereas 
if we take 

D = {u = p-'^v : V G C°^(U),p^-"nijdfC''u = on dU for k = l,2,.. .}, 

we have the self-adjoint extension corresponding to Neumann boundary conditions. Finally, 
taking 

D = {u = p-'^v-+p'^v+ : v± G C°^(U), {p^-''nijdf+p''^)/:''u = on dU for = 0, 1, . . .}, 

gives the self-adjoint extension corresponding to the Robin boundary conditions. It is 
straightforward to check that C is positive and symmetric in all three cases, provided for 
the Robin case we take /3 > 0. 

As a consequence, we can apply the functional analytic machinery of essentially self- 
adjoint operators to C in these situations. In fact, we can do better than this, based 
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on the close analogy with the finite case. It can be shown that £ with homogeneous 
Dirichlet, Neumann or Robin boundary conditions has a countable set of eigenvalues with 
corresponding eigenfunctions, smooth in the interior of U , which form an orthonormal basis 
for L^{U). This result comes from first establishing that H}{U) is compactly embedded in 
L^{U), and using this fact to apply the Fredholm alternative to a suitably chosen compact 
operator. 
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